Non-spherical proton shape and hydrogen hyperfine splitting * 
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We show that the non-spherical charge distribution of the proton manifests itself in hydrogen 
hyperfine splitting as an increase (in absolute value) of the proton Zemach radius and polarization 
contributions. 
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I. INTRODUCTION 

Hydrogen hyperfine splitting (hfs), which is predom- 
inantly due to the interaction of electron and proton 
magnetic moments, is an interesting observable. It pro- 
vides not only precision tests for quantum electrodynam- 
ics, but also valuable information on proton structure 
and strong interactions. According to Fermi's theory 
the energy difference between the two hydrogen hyperfine 
states is 



3 7T 



<% e + m v 



(1) 



where a is the fine structure constant, m e is the mass 
of the electron and fi e its magnetic moment; m p and \i v 
are the proton mass and magnetic moment. Numerically, 
this gives Ep — 1, 418, 840.11 kHz using recent values for 
these fundamental constants Q. 

The interplay between experiment and theory has been 
particularly fruitful in the case of hydrogen hfs. For ex- 
ample, when the original measurement [3J exceeded the 
prediction of Eq.JT]) by 0.26 %, which was far more than 
the experimental accuracy at that time, this stimulated 
the first quantum electrodynamics (QED) calculation of 
the electron anomalous magnetic moment [4J . The latter 
was then identified as the major reason for the ~ MHz 
discrepancy between experiment and Fermi's formula. 
Further studies revealed that hydrogen hyperfine split- 
ting is not only affected by global proton properties such 
as its mass, charge, and magnetic moment, but is also 
sensitive to the details of the spatial charge and current 
distributions in the proton and its excited states [H, Q . 
These nucleon structure effects contribute to hfs at the 
level of several tens of kHz. 

In the meantime, the transition frequency between the 
two hyperfine states has been measured [7J] with an ex- 
perimental uncertainty of about 1 mHz 



f hfs 

exp 



1, 420, 405, 751.7667 ± 0.0009 [Hz], (2) 
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corresponding to a relative accuracy of 10 -12 . Thus, hy- 
drogen ground state hfs is one of most precisely measured 
physical quantities. On the other hand, the accuracy of 
present calculations of nucleon structure effects in hydro- 
gen hfs is at best of order 10 -6 and hence many orders 
of magnitude lower than the precision with which the 
fundamental constants and QED corrections are known. 
Therefore, hydrogen hfs can be used as a high precision 
probe for investigating fine details of proton structure. 

Another way of studying the structure of the nucleon 
is electron-nucleon scattering. Employing polarized elec- 
tron beams and hydrogen targets, it has recently become 
possible to experimentally determine the p — > A + (1232) 
charge quadrupole transition form factor Q. It has 
been proposed that the quadrupole excitation of the nu- 
cleon N(939) to the A(1232) resonance is closely related 
to a quadrupole deformation of the nucleon's ground 
state charge distribution as reflected by a positive in- 
trinsic quadrupole moment Q and an intrinsic charge 
quadrupole form factor [l(| • This interpretation is based 
on relations between nucleon ground state and N — > A 
transition form factors that follow from broken strong 
interaction symmetries. 

The purpose of this paper is to explore in which way 
and to what extent the proton's non-spherical charge 
density affects hydrogen hfs. As a spin 1/2 particle the 
proton does not have a spectroscopic quadrupole moment 
and its non-spherical charge distribution does not result 
in a quadrupole interaction term in the hfs energy shift 
formula [lo| . Nevertheless, deviations from a spherically 
symmetric proton charge distribution are detectable in 
hydrogen hfs via the quadrupole excitation of the nucleon 
to the A + (1232) resonance by the atomic electron (polar- 
ization shift) and via their effect on the electromagnetic 
size of the proton (Zemach radius). 

The paper is organized as follows. In sect. [IT] we dis- 
cuss the electromagnetic N — > A transition and what we 
have learned from it about the geometric shape of the 
nucleon in some detail. In sect. IIIII the implications of 
the nucleon's non-sphericity for the hydrogen atom hy- 
perfine splitting are investigated. The paper closes with 
a summary of our results and an outlook. 
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FIG. 1: Left: Probing of nucleon structure via elastic electron- nucleon scattering eN — > e' N' involving the exchange of a 
single virtual photon 7 of four-momentum Q = (y, — q) with v and q being its energy and three-momentum transfer. Nucleon 
structure information is contained in the magnetic dipole form factor G^j(Q 2 ) and charge monopole form factor Gq{Q 2 ). 
Right: Inelastic electron-nucleon scattering eN — > e A — > e'N'ir (electro-pionproduction). The electromagnetic excitation of 
the A(1232) resonance is described by three transition form factors G^f^ A (Q 2 ), G^2* A (Q 2 ), and Gq^* a (Q 2 ). 



II. ELECTROMAGNETIC N A TRANSITION 
AND NUCLEON SHAPE 



B. Spin-flavor symmetry and electromagnetic form 
factor relations 



A. Elastic and ineleastic electron scattering 



Nucleon structure information is encoded in two elas- 
tic electromagnetic form factors, namely the charge 
monopole Gq(Q 2 ) and magnetic dipole G^(Q 2 ) form 
factors as indicated by the black dot in Fig. [1] (left). 
These form factors have been measured in elastic 
electron-proton and electron-deuteron (neutron) scatter- 
ing experiments performed at various laboratories. In 
particular, it has been shown that the proton has a fi- 
nite charge radius of about 0.9 fm In addition, the 
Fourier transforms of the elastic form factors have pro- 
vided information on the radial variation of the charge 
p(r) and current j(r) densities of the proton [l2T ]. 

Inelastic electron-proton scattering with the produc- 
tion of a single pion (electro-pionproduction) has revealed 
that the proton has a rich spectrum of excited states [l3| . 
Its lowest lying excited state with spin 3/2 and isospin 
3/2, the A(1232) resonance, plays a special role because 
it has the largest production cross section and its prop- 
erties are most closely related to those of the nucleon 
ground state iV(939). Parity invariance of the electro- 
magnetic interaction and angular momentum conserva- 
tion restrict the N — > A excitation to magnetic dipole 
(Ml), electric quadrupole (E2), and charge or Coulomb 
quadrupole (C2) transitions with corresponding transi- 
tion form factors as depicted by the large black dot in 
Fig. [1] (right). The nonzeroness of the C2 form fac- 
tor indicates that the nucleon charge distribution is not 
spherically symmetric |8[ but has an angular dependence 
p(r) = p{r, Q, $). In the following we review the connec- 
tion between the quadrupole excitation of the A(1232) 
resonance and nucleon ground state deformation using 
strong interaction symmetries as a guide. 



Aside from SU(2) isospin and SU(3) flavor symmetries, 
strong interactions are also approximately invariant un- 
der the higher SU(6) spin-flavor symmetry. The latter 
unites the spin 1/2 flavor octet baryons (2x8 states), 
among them the familiar proton and neutron, and the 
spin 3/2 flavor decuplet baryons (4x 10 states), including 
the four A(1232) states into a common 56-dimensional 
mass degenerate supermultiplet We now under- 

stand that the underlying field theory of strong inter- 
actions, quantum chromodynamics (QCD), possesses a 
spin- flavor symmetry which is exact in the large N c limit, 
where N c denotes the number of colors. Moreover, for fi- 
nite N c , spin- flavor symmetry breaking operators can be 
classified according to the powers of 1 /N c associated with 
them. This leads to a perturbative expansion scheme for 
QCD processes that works at all energy scales [151 ] . 

For N c = 3 we may just as well employ a parametriza- 
tion method (l6j . which incorporates SU(6) symmetry 
and its breaking similar to the 1/N C expansion. The ba- 
sic idea is to write for the observable under investigation 
the most general spin-flavor operator. Generally, this is 
a sum of one-, two-, and three-quark operators in spin- 
flavor space multiplied by a priori unknown constants 
which parameterize the orbital and color space matrix 
elements, and which are determined from experiment. A 
multipole expansion of the nucleon charge density oper- 
ator p in spin-flavor space up to quadrupole terms leads 
to the following invariants 

3 / 3 3 

P=P[i] +P[2] +P[3] = A^e.l- LB^ei + C^efel 



2 <Ti • CTj — (3<T izCj z — tTi ■ CTj) 
spin scalar spin tensor 



(3) 



where cti z is the z-component of the Pauli spin matrix of 



3 



c£r A (Q 2 ) 



Q 



N^A 



: = G 



N—>A 
CI 



(°) = -m 



quark i, and e, = |(1 + 3Ti Z ) is the quark charge where the elastic neutron charge form factor Gq(Q 
T iz is the third component of the Pauli isospin matrix. 
The constants A, B, and C contained in the one-, two- 
and three-quark charge density operators parametrize the 
orbital and color matrix elements so that p is only an op- 
erator in spin- flavor space. The factors 2 and -1 in front 
of the spin scalar and spin tensor terms are dictated by 
group theory and reflect the fact that both terms arise 
from a common SU(6) spin-flavor tensor [10]. In terms 
of fundamental processes (see Fig. [5]), the one-quark op- 
erator in Eq.© represents valence quarks whereas the 
two-, and three-quark operators provide an effective de- 
scription of the quark-antiquark degrees of freedom in 
the N and A. For the two-body current in Fig.[2^b) this 
is explained in more detail in Ref. (l7j . 



r C2 



Q 

1 

l r A 

10 rl 



(6) 




i Bsaaszr 



(b) 



i w aaaam 




(0 



FIG. 2: Fundamental photon-quark processes contributing to 
the form factors in Fig.fi] (a) one-quark current (p[i], j[i]), (b) 
two-quark gluon exchange current (ppjjp]), ( c ) three-quark 
gluon exchange current (p[3],j[3j). 



Evaluating the char ge o perator in Eq. ([3]) between spin- 
flavor wave functions [18| of the 56 dimensional SU(6) 
ground state multiplct, in particular, for the neutron and 
between the initial proton and final A + states, leads to 



which is valid for both the p — > A + and n — > A transi- 
tions, and that the quadrupole transition radius is deter- 
mined by the fourth and second radial moments of the 
neutron charge distribution. 

In addition, SU(6) spin-flavor symmetry leads to the 
following relations [1J| between the neutron ground state 
and the N — ► A magnetic form factors G M ~J* A (Q 2 ) = 
— v2 G M (Q 2 ), and at Q 2 = between the neutron and 
transition magnetic moments hn^a = — V^Mn- With 
the help of Eq.© and the magnetic form factor relations, 
the C2/M1 ratio in electromagnetic A(1232) excitation 
can be expressed in terms of the neutron elastic form 
factors as follows I2l| 



q\m N G^ A (Q 2 ) _ [qjmg G n c {Q 2 ) 
2Q 2 G n M (Q 2 y 
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fN—>A 
"Ml 



(Q 2 ) 

2 



(to a - m N ) m N r\ 
12 /J, 



(7) 



where |q| is the modulus of the photon's three- 
momentum and where mjv, tua are the nucleon and A 
masses. The theoretical uncertainty of this relation is 
mainly due to third order SU(6) symmetry breaking op- 
erators of order l/N 2 (three-quark currents) violating the 
magnetic form factor relation. It is estimated that such 
correction terms could lead to 10% decrease of \C2/M\\. 



rl = (56 n |p [2] +P[ 3 ]|56„ 
1 P ^A+ = (56 A +|P[2] + P[3] 156 



= 4(B-2C), 
,„ = 2V2(B- 2C).(4) 



Note that in this approach, valence quarks make no con- 
tribution and both observables are governed by quark- 
antiquark degrees of freedom in the nucleon. Hence, the 
following relation between the transition quadrupole mo- 



ment Q p ^a+ 
tained 



and the neutron charge radius r„ is ob- 



C. Comparison with experiment 

In Fig. [31 the experimental C2/M1 ratio as measured 
in electro-pionproduction is shown. The full curve rep- 
resents the Maid 2007 analysis [13] of the world C2/M1 
data, while the dashed-dotted line is based on Eq.([7]) 
which relates the inelastic N — > A and the elastic neu- 
tron form factors. For definiteness we use a Galster 
parametrization [231 ] 



V2' 



(5) 



This relation was originally derived in the constituent 
quark model with two-quark exchange currents [lj| and 
shown to hold after including three-quark operators [2(| • 
It was found that Eq.© is the zero momentum transfer 
limit of a more general relation [2l[ between the N — > A 
charge quadrupole transition form factor Gq^* a (Q 2 ) and 



g c (Q 2 ) 

G n M (Q 2 ) 



1 + dT 



GIAQ 2 



Q 2 



A AI 



(8) 



for the experimental neutron charge form factor, where 
Gm(Q 2 ) is a dipole representation of the neutron mag- 
netic form factor with Am being the dipole mass. The pa- 
rameters a and d are related to the second and fourth mo- 
ment of the neutron charge distribution respectively (2lT ] 
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and have the numerical values a = 0.9 and d = 1.75. As 
is clear from Fig. [3] (left), our theory agrees well with the 
C2/M1 data at low momentum transfers. In particular, 
at the real photon point Q 2 — we get from Eq.([7]) us- 
ing the experimental neutron charge radius and magnetic 
moment C2/M1 = —0.035, which is in good agreement 
with determinations of this ratio based on the experimen- 
tal E2/M1 ratio measured in photo-pionproduction and 
Siegert's theorem relating E2 and C2 form factors 

Moreover, the transition quadrupole moment has been 
extracted from photo- and electro-pionproduction data as 
Qp_ +A +(exp) = -0.108(9) fm 2 ^ and Q p ^ A +(exp) = 
—0.0846(33) fm 2 (25|, in good agreement with the pre- 
diction Q p ^ A + (theory) = -0.0820(20) fm 2 based on 
Eq. ([5]) and the experimental neutron charge radius r 2 = 
—0.1161(22) fm 2 . Concerning the transition quadrupole 
radius, we obtain from Eq.® and experimental values for 
the radial moments [21| of the neutron charge distribu- 
tion rc2 = 1-43 fm. The approximate equality rc2 ~ f w , 
where r w is the pion Compton wavelength suggests that 
tq2 measures the spatial extension of the qq pair distribu- 
tion in the nucleon. It would be interesting to determine 
this radius experimentally. 

At higher momentum transfers shown in Fig. [3] (right) 
the extraction of individual electromagnetic multipoles 
from the raw cross section data is more difficult as is ev- 
ident from the difference between the Jefferson Lab [26[ 
and Maid 2007 [22j analyses of the same raw data. The 
latter analysis, which leads to smaller |C2/M1| (filled cir- 
cles) than the former (open triangles) is in much better 
agreement with our theory (dashed-dotted curve). Ex- 
trapolating Eq.© to Q 2 -> oo we find [2l[ 

C2/Afl(g 2 ^oo) = -i^(^). (9) 
4 mi V a ) 

consistent with perturbative QCD, which states that 
C2 /Ml asymptotically approaches a small negative con- 
stant. 

D. Intrinsic quadrupole form factor of the nucleon 

For an interpretation of these results it is important 
to distinguish between the spectroscopic and intrinsic 
quadrupole moment of a particle. It is known that a van- 
ishing spectroscopic quadrupole moment does not neces- 
sarily imply a spherically symmetric charge distribution. 
For deformed spin and spin 1/2 nuclei this has led to 
the general concept of an intrinsic quadrupole moment, 
which can be defined for different nuclear models. The 
notion of an intrinsic quadrupole moment allows to inter- 
pret measurable quantities such as transition quadrupole 
moments in terms of the geometric shape of the ground 
state. 

The geometric shape of a spatially extended particle is 
determined by its intrinsic quadrupole moment, 

Qo= fd 3 rp(r)(3z 2 ~r 2 ), (10) 



which is defined with respect to the body-fixed frame. If 
the charge density is concentrated along the z-direction 
(symmetry axis of the particle), the term proportional to 
3z 2 dominates, Qq is positive, and the particle is prolate 
(cigar-shaped). If the charge density is concentrated in 
the equatorial plane perpendicular to z, the term propor- 
tional to r 2 prevails, Qq is negative, and the particle is 
oblate (pancake-shaped) . 

We calculated the intrinsic quadrupole moment of 
the proton and A + in the quark model including two- 
body exchange currents Q , which effectively describe the 
quark-antiquark degrees of freedom in the nucleon, and 
found 

Ql = ~rl Qf=rl (11) 

Thus, the intrinsic quadrupole moment of the proton is 
given by the negative of the neutron charge radius and is 
therefore positive, whereas the intrinsic quadrupole mo- 
ment of the A + is negative. This corresponds to a prolate 
proton and an oblate A + shape. The model results also 
suggest that the nonsphericity of the proton charge den- 
sity is mainly connected with collective quark-antiquark 
degrees of freedom, the distribution of which has a pro- 
late shape. 

The concept of an intrinsic quadrupole moment of the 
nucleon can be generalized to an intrinsic quadrupole 
charge distribution and a corresponding form factor [l0( . 
To show this, we first decompose the proton and neu- 
tron charge form factors in two terms G sym and Gdef, 
coming from the spherically symmetric and the intrin- 
sic quadrupole part of the physical charge density repec- 
tively 

G P C (Q 2 ) = G% m {Q 2 )- l -Q 2 G def {Q 2 ), 
GUQ 2 ) = G n sym (Q 2 ) + l -Q 2 G def {Q 2 ). (12) 

The factor Q 2 in front of Gdef arises for dimensional rea- 
sons and guarantees that the normalization of the charge 
form factors is preserved. In coordinate space this cor- 
responds to the usual multipole decomposition of the 
charge density 

p(r) = po(r)r o °(r) + p 2 (r)Y 2 (v) + (13) 

monopolc quadrupole 

where the po(r) part gives rise to G sym (Q 2 ) and the p2{r) 
part is connected with Gd e f(Q 2 )- In terms of fundamen- 
tal photon-quark processes Gdef arises mainly from two- 
and three-quark currents. 

For the intrinsic charge quadrupole form factor 
Gdef{Q 2 ) we find 

Gdef(Q 2 ) = -V2G^ A (Q 2 ) = -^G n c (Q 2 ), 
Gde f (0) = -rl = Ql (14) 
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FIG. 3: The C2/M1(Q 2 ) = Si+/M 1+ (Q 2 ) ratio for low (left) and high (right) four-momentum transfers. The full curve is a 
fit of the experimental C2/M1 ratio as determined from the world electro-pionproduction data. The dashed-dotted curve is 
calculated using the form factor relation of Eq.(0. The open triangles and the dashed curve are from a previous analysis of 
the same data [2fj]. Figure taken from Ref. [22 ]. 



where we have used Eq.([6]). The zero momentum limit 
follows from 1' Hospital's rule and Eq. fTTj) . This shows 
that Gdef(Q 2 ) as denned in Eq. (fl"4|) is the proper gener- 
alization of the intrinsic quadrupole moment Qq to finite 
momentum transfers. 

To exhibit the effect of the intrinsic quadrupole form 
factor on the elastic nuclcon form factors we insert 
Eq.lfHf into Eq.([l2]) and obtain 

G p c (Q 2 ) = G% m (Q 2 )-G^Q 2 ) = G^m-G^m, 

spherical deformed 

G n c {Q 2 ) = ^ Q 2 Gdef{Q 2 ), (15) 

where the isoscalar nucleon charge form factor is defined 
as G J C S (Q 2 ) = G P C (Q 2 ) + G£(Q 2 ) and G% m = 0. Thus, 
the relation between the N — > A and neutron charge form 
factors of section III Bl is seen here to have an important 
implication for the nucleon itself, which can be summa- 
rized as: The neutron charge form factor is an observable 
manifestation and quantitative measure of the nucleon's 
intrinsic quadrupole form factor. The latter manifests 
itself also in the proton charge form factor. 

There are several observable consequences of Eq. fTl)) 
and Eq. (fl~5|) as discussed in Ref. [l(|. At low Q 2 the 
nucleons's prolate deformation is reflected in a proton 
charge radius increase by an amount — r„, or more di- 
rectly by a newly introduced size parameter r^ e y = r^, 2 = 
(7/10) (r^/r 2 ) that can be experimentally determined. 
At intermediate Q 2 it leads to the conclusion that the dip 
structure observed in the proton charge form factor [12j | 



at around Q 2 ~ 0.2 GcV 2 is due to a corresponding struc- 
ture in the neutron charge form factor at the same Q 2 . 
Finally, at high Q 2 it explains the observed decrease of 
the charge over magnetic form factor ratio (27| . 

We close this section by stating that our introduction 
of an intrinsic quadrupole moment and quadrupole form 
factor of the nucleon should be viewed as an attempt to 
explore the consequences of the experimental sign and 
size of the N — > A quadrupole transition form factor 
for nucleon ground state structure, which in turn has a 
bearing on hydrogen hyperfine splitting. 



III. HYDROGEN HYPERFINE SPLITTING 
AND NON-SPHERICAL PROTON SHAPE 

A. Fermi energy 

It is well known that the ground state energy of atomic 
hydrogen is split into two levels due to the interaction 
of the magnetic moments of the electron fj, e and proton 
fi p , which can be either aligned (F = 0) or antialigncd 
(F = 1). The energy difference between these two states 
is of order 10~ 6 eV which is small compared to the -13.6 
eV binding energy of the ground state. For spherically 
symmetric electronic S states, the magnetic dipole-dipole 
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interaction Hamiltonian can be expressed as 
2 

H F = --fi p - fi e S {3) (r p - r e ) 

where r p and r e are the proton and electron position 
coordinates and 5^ is the three-dimensional Dirac 5 
function. The second equation follows after rewriting 
the magnetic moments in terms of spin operators as 
H P = (1 + k) [e/ (2m p )] (T p and /x e = (1 + a) [e/ (2m e )] <x e . 
By convention the proton anomalous magnetic moment 
1 + k is included in the Fermi energy, while the anoma- 
lous electron magnetic moment 1 + a is included in the 
QED corrections introduced in sect. IIII 51 

After taking matrix elements of Eq. (fH)]) between hy- 
drogen ground state electron wave functions v E' e (r e ) one 
obtains for the hyperfine level splitting 

-(^■(T e ) F=0 ^ e (0)| 2 , (17) 

where (5' 3 '(r p — r e ) has been evaluated for r e = r p = 
(point nucleon limit). We use the standard normalization 
N of the hydrogen atom ground state wave function 

\* e (0)\ 2 = [4=42 J = N *> «fl = — » (18) 

where as is the Bohr radius, which is given in terms 
of the reduced mass m r — m p m e /(m p + m e ) of the 
electron-proton system and the fine structure constant 
a = e 2 /(47r). The spin matrix element in Eq. (|17p gives 
< cr p ■ cr e >f—\= 1 and < cr p ■ cr e >f=o= — 3 so that 
one obtains the Fermi energy formula of Eq. (Q} and the 
numerical value E F = 1,418,840.11 kHz. 



B. QED and proton structure corrections to the 
Fermi energy 



to be compared with the experimental value in Eq.([2]). 
One readily notices that there is still a discrepancy be- 
tween theory and experiment, namely Eq^q — E^ xp s = 
46.29 kHz or 32.63 ppm. 

Second, adding nucleon structure dependent relativis- 
tic recoil corrections S rec — 5.85 ppm to the theory [3(| 
this discrepancy increases to 38.48 ppm, i.e., a significant 
deviation by which the theoretical value exceeds the mea- 
sured one. Third, as a consequence of the proton's finite 
size, its magnetic moment is distributed over an extended 
spatial region. This weakens the magnetic interaction 
with the atomic electron and reduces the hfs. A first es- 
timate of the proton size effect can be obtained from an 
expansion of the electron wave function for small radial 
distances v f , e(^) = -/Vexp(— r/oe) = N (1 — r/os + • • ■ )■ 
The spatial extension of the proton magnetic moment 
distribution is then taken into account by evaluating the 
electron wave function for r e = r m ^ 0, where r m is the 
proton magnetic radius. This provides a correction term 
to the Fermi energy of the form [3l| 

Eprofonsize = E F (l ~ 2 (20) 

and the numerical estimate of the proton size effect 
-2r m /a B ~ -2 ■ 10~ 5 A/0.5l = -40 ppm. When these 
corrections are added to Eq. (fP9|) . one obtains a reduc- 
tion of the theoretical hfs which is of the right size to 
achieve agreement between theory and the experimental 
value in Eq.((2|) at the ppm level. The recoil and finite 
size corrections originate both from second order elastic 
electron-nucleon scattering depicted by the two-photon 
exchange diagram in Fig. [4] (left) and are conventionally 
denoted as S rec and Sz- 

Fourth, there are also inelastic contributions (nucleon 
polarization), which involve intermediate excited pro- 
ton states, e.g., the A(1232) resonance as shown in 
Fig. d] (right). The elastic and inelastic second order 
electron-proton interaction terms are generically referred 
to as proton structure contributions. The most impor- 
tant corrections to the Fermi energy are then 

E theory = E F (1 + &QED + Srec + $Z + Spot) ■ (21) 



It has become customary to express QED and proton 
structure corrections as parts per million of the Fermi 
energy, i.e., 1 ppm =10~ 6 i?F=l-41884 kHz. The most 
important corrections to the pointlikc dipole-dipolc inter- 
action energy in Eq.|T]) are due to (i) QED, (ii) nucleon 
recoil, (iii) finite proton size, and (iv) proton polarization 
effects [28]. First, there is the anomalous magnetic mo- 
ment of the electron, which is mainly caused by the QED 
vertex correction. This and other smaller QED contribu- 
tions lead to a 1136.1 ppm increase of the theoretical 
hfs 

EqEE = E f (1 + 8 QE d) = 1,420,452.04 kHz (19) 



C. Zemach radius and proton shape 

A careful analysis of the proton finite size correc- 
tion in hydrogen hyperfine splitting was performed by 
Zemach [5j. Assuming rigid (unpolarizable) , spherically 
symmetric charge and magnetization distributions for the 
proton, the following expressions were derived 

n Jo Q L Mp 

Sz = -2r z /a B (22) 



7 




FIG. 4: Two-photon exchange diagrams (crossed photon diagrams are not shown) from which the elastic (Zemach and recoil) 
and inelastic (polarization) nucleon structure corrections to hydrogen hyperfine splitting are derived. The large black circles 
represent elastic (left) and inelastic (right) nucleon electromagnetic form factors as measured in electron-nucleon scattering 
depicted in Fig. Q] 



where Gq(Q 2 ) and G V M (Q 2 ) are the elastic charge and 
magnetic form factors of the proton. In contrast to 
the estimate in Eq.(|20|) the Zemach correction S z de- 
pends on the details of both, charge and magnetic mo- 
ment distributions. The term —1 in the integrand of 
Eq. (f2"2")l is necessary because the point nucleon limit, 
Gq(0) — G p M (0)/fi p — 1, is already included in the 
Fermi energy Ep and must be subtracted to avoid double 
counting. There is also a radiative correction term S ra d 
due to electronic vacuum polarization [32| which changes 

Sz — * Sz (1 + S ra d)- 

In order to separate the effect of proton's non-spherical 
charge distribution on the Zemach radius we insert for 
G P C (Q 2 ) the decomposition of Eq.([T5|) into Eq.([22"f and 
obtain 



rz 



'dQ 



jlp 



spherical 



deformed 



rz(sym) + r z {def). 



(23) 



Thus, the proton Zemach radius is decomposed into two 
terms coming from the spherically symmetric and non- 
spherical parts of the proton charge distribution respec- 
tively. 

For an estimate of the spherically symmetric contribu- 
tion we assume dipole forms for the isoscalar charge and 
proton magnetic form factors and obtain the compact 



two-parameter formula 



r z (sym) 



3 A) + 9 A? A M + 11 A? A 2 M + 9 A/ A\ 



M 



3A 4 M 



A/Am (A/ + A M ) 3 



which hitherto has not appeared in the literature. Here, 
A/ and Am are the inverse size parameters of the isoscalar 
charge and proton magnetic form factors. They are re- 
lated to the corresponding mean square isoscalar charge 
and proton magnetic radii as A 2 = 12/rf with rf = 



r c(p) + r c( n ) an d A 2 M = 12 /r 2 M (p). In the limit, 
A/ = A M = A, Eq.(|2i]) reduces to r z (sym) = 35/(8 A), 
i.e., a standard expression that has been used by sev- 
eral authors [12, [H, [38[. The contribution of the non- 
spherical part of the proton charge distribution can also 
be analytically calculated if one uses the two-parameter 
Galster form of Eq.© for G£(Q 2 ) 



rz(def) = --rl(A M m) 



16A M + 29A 2 M m + 20 A M m 2 
8 (A M + m) 4 



-,(25) 



where we have used the same dipole mass Am for the 
proton and neutron magnetic form factors. The quanti- 
ties r 2 and m 2 are related to the Galster parameters a 
and d as r 2 n — 3a/i n /(2m^) and m 2 :— Am 2 N /d. Eq. 



and Eq. (|25[) are two main results of this paper. 

Inserting experimental values for the dipole parameters 
Am and Aj, as well as for r 2 and the Galster parameter 



d = 1.75, one finds the numerical result 



rz = r z (sym) + r z (def) 
= 1.1091 fm 



1.0627 fm + 0.0464 fm 

(26) 
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corresponding to Sz = —(40.14 + 1.75) ppm = —41.89 
ppm. The proton deformation contribution is numeri- 
cally of the same size but of opposite sign as the neu- 
tron Zcmach radius [34[. The sign change is obvi- 
ous from the definition in Eq. (|23p and the approximate 
equality in magnitude rz(def) ~ \rz{n)\ follows from 
the near equality of the normalized proton and neutron 
magnetic form factors in the relevant moment transfer 
range. Including the radiative correction gives our final 
result for the proton Zemach radius 8z = —41.89 ppm 
(1+0.0151) = -42.52 ppm. With the Zemach radius con- 
tribution included one finds that the discrepancy between 
theory and experiment reduces to (38.48 — 42.52) ppm 
4.04 ppm. 

To obtain a better estimate of the effect of the proton's 
non-spherical shape on rz one would have to determine 
the proton and neutron charge form factors more accu- 
rately in the low momentum transfer region where possi- 
ble deviations from the smooth dipole and Galster fitting 
curves may affect rz(def) significantly. Electron- proton 
scattering experiments dedicated to explore this low mo- 
mentum transfer region with higher precision are being 
planned [35] ]. 



the functions Pi(rj) and /?2 (''?) are defined as 

= ^(-3r, + 2^ + 2(2-^)^ + 1)), 

ftfa) = 1 + 277 - 2^ + 1), (28) 

with r\ := v 2 /Q 2 . Furthermore, F^iQ 2 ) is the Pauli form 
factor of the proton, which is defined in terms of the 
charge monopole and magnetic dipole form factors as 
Fi(Q 2 ) = (GitiQ 2 ) - Gl (Q 2 ))/(l + r) with Ff (0) = n 
being the anomalous part of the proton magnetic mo- 
ment. Note that both terms in the integrand of 8\ di- 
verge for Q 2 = but the singularity coming from the 
second term is cancelled by an analogous singularity of 
the first term according to the Drell-Hearn, Gerasimov 
sum rule [36| . 

The polarization contribution has been calculated by 
several authors [3(| H3|- However, only the contribution 
of the magnetic dipole (Ml) transition to the A(1232) 
has been studied in some detail. There has been no 
prior attempt to calculate the contribution of the charge 
quadrupole (C2) transition to the A(1232). In order to 
investigate the effect of the latter on the polarization 
shift, we express the spin structure functions in terms 
of virtual photon absorption cross sections [301 ] 



D. Proton quadrupole polarization shift 

The proton polarization shift 6 po i in hydrogen hfs is 
caused by two-photon exchange diagrams with nucleon 
resonances as intermediate states, of which the lowest- 
lying A(1232) as shown in Fig. [4] (right) is expected to be 
most important. Here, we focus on those diagrams where 
one of the photons is a longitudinal charge quadrupole 
(C2) photon that probes the non-spherical charge distri- 
bution in both the ground and excited states. 

The polarization shift is usually defined in terms of 
integrals over the two spin-dependent structure func- 
tions 51(24 Q 2 ) and g2(x 1 Q 2 ) of the proton, where x = 
Q 2 /(2m p v) is the Bjorken scaling variable and v is the 
energy transfer carried by the virtual photon (see Fig.[T]). 
The following formulae for proton polarization in hydro- 
gen hfs have been established [3(| 



Spol 

Si 



a m e 



2tt m p (l + k) 
9 f°° dQ 2 

o m 2 r x th 



(Si + 8 2 



Q 2 Jo 



dxpx{r]) gi{x,Q ) 



Urn 2 I 

p Jo Q 4 



°°rin 2 r r x th 



dx/3 2 (r)) g 2 (x,Q ) 



,(27) 



where x t h = Q 2 / \2m p + m 2 + Q 2 ) is the threshold for 
one-pion production with m, being the pion mass, and 



91 (y, Q 2 ) 



2 _ m p u(l - Q 2 /(2m p v)) 



&Tr 2 a(l + Q 2 /is 2 ) 



°i/2(v, Q 2 ) - *3/ 2 (f, Q 2 ) - 2 - o' LT (v, Q 2 ) 



2 _ m p v(l- Q 2 /(2m p v)) 
877^0 (1 + Q z lv z ) 

(-cri/2(v, Q 2 ) + 0-3/2(1/, Q 2 ) - 2 ^ a' LT {v, Q 2 )^j .(29) 

Here, CT1/2 and 03/2 are the transverse cross sections 
with photon-nucleon helicity 1/2 and 3/2, and a' LT is 
the longitudinal-transverse interference cross section. 

For excitation energies below 1 GeV, which are most 
important for hydrogen hfs, the cross section is domi- 
nated by multipole transitions to specific nucleon reso- 
nances and can be written in terms of the helicity ampli- 
tudes A l / 2 , ^3/2, an d S'1/2 as (ioj 



01/2 
03/2 



OiT 



2-i^i/ 2 | 2 , 
2-9^3/ 2 | 2 , 



= -V2tt 



m p Q 
W\q\ 



bS*, 2 Ai/2 



(30) 



where W is the total center of mass energy, q the jN 
center of mass three-momentum and b the resonance line 
shape, which at resonance reduces to b — 2/(7rr T 7\r) with 
r^jv being the pion decay width of the resonance. In par- 
ticular, for A(1232) excitation, the helicity amplitudes 
can be expressed via the inelastic N — > A transition form 
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factors [HI shown in Fig. Q] (right) 



A 1/2 (Q 2 ) 



A 3/2 (Q 2 ) 



S 1/2 (Q 2 ) 



2m 



P v Kw 



rtN—>A(n2 

"Ml 



3tt 



2m„ V Kw 



G 



N—>A 
Ml 



(0 2 



^£2 W ) 



l\ m p s~tN—>£± li^fi 

KJ 



. 2tt q 2 „jv->a (/ ,2 

' w avT Gc2 



3V6 



(<2 ) 



(31) 



where K\v = (m 2 ^ — m N ) / {2m^) is the energy transfer 
in the center of mass frame at Q 2 = 0. 

The contribution of the N — > A charge quadrupole 
(C2) transition form factor to the polarization shift, 
which comes solely from the a' LT part in Eq. ([29|) can 
now be evaluated. We obtain with m^ = m p 



S^ A (C2) = 

2 mr 



A 



~'dQ (v\ 



GUQ 2 ) 



G n M {Q 2 ) 



5^ A (C2) 



ft 77l N u 

ml 



/?2 



Q 2 J Q 



,(32) 



-Q 2 ) / (2mAr) is the resonance en- 



where = {m\ — m 2 N 
ergy in the laboratory frame. To derive Eq. (j32|) we have 
made use of the form factor relations of sect. Ill Bl and 
neglected the small E2 x G2 contribution. Numerically, 
we obtain with values for Am and m as in sect. MI Cl and 
r^jv = 0.12 GeV the following estimates 



<5f^ A (<72) = 
5^ A {C2) = 

Cr A (C2) = 



-0.994 

+0.292 
—0.16 ppm. 



(33) 



The negative sign of the quadrupole polarization 
shift is indicative of a prolate (cigar-shaped) intrinsic 
quadrupole deformation of the proton's charge distribu- 
tion. Interestingly, the charge quadrupole (C2) polariza- 
tion contribution in Eq. ([33|) is not small compared to the 
magnetic dipole (Ml) polarization 5™ J* A (M1) = -0.12 
ppm derived earlier (37l. [38^. Therefore, it will be impor- 
tant to also calculate the polarization shift induced by 



transverse electric (E2) excitation of the A(1232). Fi- 
nally, we also give our result for the Ml contribution 
coming from the convergent g 2 integral b~^~ * A (_M1) = 
—0.69 ppm compared to 8 2 = —0.42 ppm including ad- 
ditional resonances [371 ] . 

IV. SUMMARY 

The nonzeroness of the empirical p — > A + (1232) 
quadrupole transition form factor provides evidence that 
the charge distribution of the proton ground state devi- 
ates from spherical symmetry. Employing SU(6) spin- 
flavor symmetry as a guide we have derived a relation 
between the p — ^ A + quadrupole transition and neutron 
charge form factors. It has been shown that this relation 
agrees with the experimental data from low to high mo- 
mentum transfers. On this basis, we have proposed that 
the proton can be assigned an intrinsic charge quadrupole 
form factor and a positive intrinsic quadrupole moment 
corresponding to a prolate (cigar-shaped) distribution of 
the proton charge. 

We have then investigated how hydrogen ground state 
hfs is affected by the proton's non-spherical charge dis- 
tribution. We have shown that the latter is reflected in a 
positive deformation contribution to the proton Zemach 
radius, where the increment is given by the modulus of 
the neutron Zemach radius. A second consequence of 
the proton's prolate shape is a polarization shift due to 
the p — ► A+(1232) charge quadrupole transition that is 
quantitatively described by the neutron charge form fac- 
tor. We have presented a numerical estimate for this term 
and found that it provides a negative contribution to the 
polarization shift in atomic hydrogen hfs that exceeds the 
one coming from the previously calculated p — > A + (1232) 
magnetic dipole transition. 

In view of these results it would be interesting to also 
explicitly calculate the polarization shift induced by the 
transverse electric quadrupole form factor G I ^^* A . Be- 
cause the experimental accuracy of hfs measurements ex- 
ceeds the theoretical accuracy by several orders of magni- 
tude, hydrogen hfs will remain a high precision probe for 
proton structure for many years to come. It is quite pos- 
sible that the interplay between hfs experiment and the- 
ory will provide independent evidence for the quadrupole 
deformation of the proton's charge distribution. Con- 
ceivably, hydrogen hfs is also sensitive to higher mag- 
netic multipoles, in particular to an intrinsic magnetic 
octupole term in the proton's spatial current distribu- 
tion. We hope to discuss these matters in a future com- 
munication. 
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